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摘 要 :为 了 提高 聚 醚 栈 酮 (polyether ether ketone , PEEK ) 材料 颈椎 椎 间 融 合 器 的 融合 率 , 本 研究 利 
用 有 限 元 方法 对 其 进行 三 维 结构 的 拓扑 优化 设计 。 首 先 建 立 C3 ~ C7 节 段 的 颈椎 有 限 元 模型 并 验 
证 ,然后 创建 C5 ~ C6 节 段 椎 间 融 合 器 植 入 和 前 路 钢板 内 固定 的 模型 。 将 PEEK 融合 器 设置 为 椭 
\O 〇 盘 型 和 箱 型 两 种 初始 设计 域 。 先 后 进行 单 工 况 和 多 工 况 下 融合 器 拓扑 优化 ,并 结合 加 权 柔 度 优化 。 

月 经 过 两 次 拓扑 优化 后 ,再 对 融合 器 的 构 型 提取 分 析 。 研 究 结 果 表 明 :PEEK 融合 器 的 优化 结构 为 环 

刁 形 ,前 方 缺 失 , 材 料 集中 分 布 于 后 部 及 两 侧 ; 当 体积 分 数 分 别 为 50% 和 30% 时 ,融合 器 优化 为 多 孔 
一 径 和 单一 大 孔径 的 形状 ,前 者 比 后 者 增加 了 与 终 板 的 接触 面积 ;优化 后 的 融合 器 通过 降低 自身 应 
2 力 ,增加 了 内 部 植 骨 块 所 受 的 应 力 。 本 研究 通过 加 权 柔 度 的 方法 构建 的 颈椎 融合 器 优化 模型 ,不 受 
CD 初始 设计 域 变 化 的 影响 ， 且 可 以 有 效 减少 应 力 遮挡 ,增加 植 骨 块 的 应 力 传导 ,将 有 助 于 促进 融合 。 


: 忆 关 键 词 :颈椎 ;融合 器 ;拓扑 优化 ;有 限 元 分 析 ; 聚 醚 栈 酮 
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Abstract :The three-dimensional topological optimization of cervical interbody fusion cage was developed 
by the finite element method to enhance the fusion rate of the cervical interbody cage. The finite element 
(FE) model of C3-C7 segments was constructed and validated. The fusion FE model was developed by in- 
serting the cage into the C5-C6 disc space and internal fixation used anterior plate-screw system. The ini- 


tial structures of cages made of polyether ether ketone ( PEEK) were set up into two types of elliptic cage 
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and box. Topological optimization and weighted compliance optimization were developed under single and 


multiple loading conditions. The structural configuration analysis of the cage was carried out after two times 


of topological optimization. The results showed that the optimal structure of the cage is annular, and the 


materials of the cage are concentrated in the back and both sides;the cage is optimized for multi-pore 


shape and single large-pore shape when the volume fraction is set to 5096. and 3046 , respectively , and the 


former has more contact area with the endplate than the latter;the optimized cage increases the stress on 


the internal bone graft by reducing its own stress under loading. In this study , the optimal PEEK cage was 


designed by topological and weighted compliance optimization method, which was not affected by the 


change of initial structures of cage. It can effectively reduce the stress shielding and increase the stress 


stimulation of internal bone graft, which will help to promote the cervical fusion. 
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由 于 经 济 发 展 导致 的 工作 压力 增 大 ,颈椎 病 的 
发 病 率 逐 年 升 高 ,严重 的 病例 可 造成 神经 根 和 /或 疹 
BE IE, 因而 需 接受 手术 治疗 。1955 年 , Uschold 
等 四 开创 了 颈椎 前 路 间 盘 切除 内 固定 术 (anterior 
cérvical discectomy and fixation, ACDF)。 目 前 ,该 术 
式 也 成 为 治疗 各 种 类 型 退行 性 颈椎 疾病 的 最 常用 手 
术 汶 法 。 近 年 来 , 椎 间 融 合 器 的 使 用 已 成 为 ACDF 
季 染 的 重要 组 成 部 分 , 它 有 助 于 实现 术 后 的 即刻 稳 
污 3 振 开 并 维持 椎 间 孙 高度, 已 经 取得 了 满意 的 临床 
疗效 2” 。 尽 管 如 此 , 随 着 临床 手术 病例 的 增多 ,与 融 
你 器 直接 相关 的 并 发 症 也 时 常 出 现 ,最 常见 的 就 是 
界面 间 融 合 不 良和 融合 器 下 沉 ” 。 前 者 会 造成 融合 
器 移 位 ， 而 后 者 则 会 引起 手术 节 段 的 椎 间 队 高度 丢 
灾 写 两 者 均 有 导致 湖 髓 和 /或 神经 根 受 压迫 的 风险 ， 
从 是 影响 ACDF 术 后 患者 的 疗效 ， 其 至 导致 手术 失 


并 通过 加 权 柔 度 的 方法 构建 颈椎 融合 器 的 优化 模 
型 最 后 对 优化 后 的 融合 器 构 型 进行 生物 力学 性 能 


(m. 


1 材料 与 方法 


1.1 颈椎 三 维 有 限 元 模型 的 建立 与 验证 


1.1.1 实验 对 象 

本 研究 选取 健康 成 年 男性 ,年龄 32 岁 , 身高 
170 cm ,体质 量 67 kg, 既往 无 脊柱 间 盘 退 变 、 畸 形 、 
骨折 及 肿瘤 病史 , 且 骨 密度 正常 。 患 者 取 仰 卧 位 , 采 
用 64 排 螺 旋 CT( Brilliance 64 Philips) 扫描 C3 ~ C7 
节 段 ,设置 像素 大 小 为 0.328 mm x0. 328 mm 图 像 
厚度 为 0.625 mm ,将 图 片 以 DICOM 格式 导出 。 


败 公 引发 医疗 纠纷 。 因 此 ， 如 何 通过 对 颈椎 融合 器 
的 多 化 设计 来 降低 此 类 并 发 症 的 发 生 率 ， 已 成 为 背 
柱 外 科 领 域 迫 切 需要 解决 的 问题 。 尽 管 聚 醚 醚 酮 
(polyether ether ketone, PEEK) 材料 以 其 优良 的 生 
物力 学 特性 ， 已 取代 铁合金 而 被 广泛 应 用 于 疹 柱 融 
合 器 的 制备 ， 但 与 相 邻 椎 体 仍 存在 生物 力学 适 配 性 
问题 。 现 阶段 对 PEEK 材料 颈椎 融合 器 的 生物 力学 
适 配 性 改进 多 集中 于 结构 的 优化 。 相 比 传统 的 尺寸 
和 形状 优化 ， 以 材料 分 布 为 优化 对 象 的 拓扑 优化 方 
法 不 但 适用 于 力学 领域 ， 也 开始 应 用 于 兰 柱 内 
植物 的 改进 设计 上 ， 但 针对 颈椎 融合 器 的 拓扑 优化 
还 罕有 报道 ， 目 前 的 研究 仅 适 用 于 腰椎 融合 器 "” ， 
且 优 化 对 象 的 构 型 单一 ， 难 以 评估 其 优化 方法 是 否 
具有 普 适 性 。 基 于 近年 来 对 贷 椎 有 限 元 建 模 及 相关 
生物 力学 分 析 的 开展 "”" ， 本 研究 旨 在 利用 有 限 元 
分 析 方 法 对 PEEK 颈椎 融合 器 结构 进行 拓扑 优化 ， 


1.1.2. 有 限 元 模型 的 建立 

1) 将 DICOM 格式 文件 导入 Mimics ( Mimics 
11. 01 , Materialise ) 软件 中 进行 椎 体 几 何 信 息 的 提 
取 , 再 导入 Hypermesh( Hyperworks 12.0,Altair) 软件 
中 建立 网 格 模 型 。 对 C3 ~ C7 这 5 个 椎 骨 模 型 进行 
网 格 划 分 后 ,再 对 所 划分 的 网 格 分 配 部 件 。 基 于 
Panjabi4 ' 的 研究 结果 ,本 研究 将 椎 骨 最 外 层 的 壳 
单元 设置 为 0.5 mm, 即 外 层 表 面 厚 0. 5 mm 的 二 维 
单元 为 皮质 骨 , 同 时 选取 上 下 0.5 mm 厚 的 二 维 单元 
为 终 板 ,其 余 内 部 实体 单元 为 松 质 骨 。 根 据 解剖 特 
性 , 骨 核 设 定 为 三 维 实心 单元 (图 1) 所 示 , 其 横 截 面 
积 约 占 整个 椎间盘 的 5096 '” , 髓 核 外 周 的 纤维 环 基 
质 也 为 实体 单元 ,并 附 有 4 层 网 状 单元 来 模拟 纤维 
软骨 。 模 型 中 还 建立 了 5 组 主要 的 颈椎 韧带 ,所 有 
T9] 352 H] Truss 单元 模拟 , 且 设 为 不 可 压缩 的 材料 


特性 。 


600 


1 C5-Ce 椎 间 纤 维 环 三 维 网 格 模型 
Fig.1 3D mesh model of C5-C6 intervertebral annulus fibrosus 
2) 模 型 装配 与 材料 属性 的 设置 。 
椎 骨 、 椎 间 盘 .韧带 都 已 在 Hypermesh 中 划分 完 


单元 网 格 ,整体 模型 如 图 2 所 示 。 模 型 中 设置 小 关 
节 为 有 限 滑 移 的 面 - 面 接触 ,摩擦 系数 为 0.01( 实际 
关节 面 滑动 较 大 , 摩 探 较 小 ,主要 靠 关 节 突 的 几何 结 
构 起 到 滑动 限制 作用 ) ” 。 模 型 中 皮质 骨 、 松 质 骨 、 
椎间盘 等 组 织 均 被 视 为 各 向 同性 、 均 质 的 线 弹 性 材 
料 ,主要 材料 属性 参数 为 弹性 模 量 、 泪 松 比 和 横 截 而 
积 有 所 有 材料 属性 的 赋值 参考 文献 1 18-21 ] ,各 部 分 
AREE 的 材料 参数 如 表 1 所 示 。 


:202303.003 


图 2 颈椎 C3 ~ C7 节 段 完整 有 限 元 模型 
SFig. 2 Finite element model of C3-C7 cervical vertebra 
Sé I 有 限 元 模型 的 材料 属性 与 单元 类 型 
a Tab. 1 
c 
eh 


The material properties andelement 


types of the finite-element models 


T 
元 dE 面 
53i 85 2 dad mex "8 dd 
椎 骨 
皮质 骨 S4 12000 — 0.29 厚度 0.4 mm 
松 质 骨 C3D4 450 0.29 - 
终 板 S4 1 200 0.4 ”厚度 0.4 mm 
关节 软骨 C3D8 10 0.4 = 
椎间盘 
BEA C3D8 1 0. 49 = 
纤维 环 基质 C3D8 3.4 0.4 一 
纤维 环 T3D2 450 0.3 0.01 
韧带 
Bü AI CALL) T3D2 40 0.3 6.0 
JE 2M CPLL) T3D2 30 0.3 5.0 
3E (CL) T3D2 5 0.3 46.0 
黄 韧 带 (FL) T3D2 5 0.3 5.0 
Ape] 97 ( ISL) T3D2 10 0.3 10.0 
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3 ) 模型 边界 条 件 .载荷 设置 

模型 前 处 理工 作 在 Hypermesh 中 完成 ,再 将 模 
型 文件 导入 Abaqus( Abaqus 6. 14-2 , Simulia) 软件 中 
进行 边界 条 件 、 载 荷 的 设置 。 参 考 体 外 经 典 实 
35 75 ,将 颈椎 CT 椎 骨 的 底部 完全 固定 (约束 底部 单 
元 节点 的 6 个 自由 度 ) ,在 C3 椎 体 上 方 的 1 个 参考 
点 施加 73.6N 的 轴 向 压缩 预 载 荷 来 模拟 头 部 重量 (此 
参考 点 为 施 力 点 ,与 C3 上 终 板 做 节点 耦合 ) 。 有 颈椎 的 
各 方向 运动 都 通过 在 参考 点 施加 相应 力矩 实现 。 


1.1.3 验证 


对 模型 依次 施加 1.5 Nm 的 力矩 后 ,测量 颈椎 前 
届 \ 后 伸 、 左 右 侧 弯 和 左右 轴 向 旋转 的 活动 度 (range 
of motion, ROM) 。 将 实验 数据 分 别 与 体外 尸体 实 
验 ” 汪 以 及 有 限 元 实验 ” 的 结果 做 对 比 (图 3 ~ 图 
5) ,从 图 中 可 以 看 出 ,在 各 个 方向 上 ,本 颈椎 模型 的 活 
动 度 及 变化 趋势 均 与 文献 中 的 数据 结果 有 较 好 的 一 臻 
性 。 因 此 ,本 研究 创建 的 颈椎 有 限 元 模型 是 有 效 的 。 


M m 文献 [20] 


1 加 文献 [22] 

p 文献 [21] 
m Ay 

10 
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C3XC4 CHC5 C5C6 C6-C7 
图 3 有 颈椎 C3 ~ C07 节 段 届 伸 关节 活动 度 
与 文献 实验 数据 对 比 结 


Fig.3 The range of motions of flexion and extension 


活动 度 /(9) 


at C3-C7 cervical vertebra were compared 


with experimental data in literatures 


活动 


C3-C4 C4-C5 C5-C6 C6-C7 
图 4 颈椎 C3 ~ C7 节 段 双 侧 弯曲 活动 度 
与 文献 实验 数据 对 比 结果 


Fig.4 The range of motions of the bilateral flexural 


at C3-C7 cervical vertebra were compared 


with experimental data in the literature 


14 m 文献 [20] 
m 文献 [22] 


CER — C4-C5 — C5-C6 —— C6-C7 
图 5 有 颈椎 C3 ~ C7 节 段 左右 扭转 活动 度 
与 文献 实验 数据 对 比 结 


Fig.5 The range of motions of the left and right torsion 
at C3-C7 cervical vertebra was compared 


with experimental data in the literature 
1.2. 颈椎 融合 器 结构 的 拓扑 优化 设计 方法 


134 35 4E 3-38 PX ACDF 模型 的 建立 

〇 ) 在 前 文 所 建立 的 C3 ~ C7 节 段 有 限 元 模型 的 基 
础 旦 ,选择 C5 ~ C6 节 段 构建 ACDF 的 手术 模型 , 因 
法 该 节 有 段 是 颈椎 疾患 发 病 与 手术 的 最 常见 节 段 *。 
杰 研 究 采用 4 枚 单 皮质 螺钉 和 锁定 钢板 , 均 为 钛 合 
< TECE 2110 000 MPa,u 20.33) ,钢板 采用 CSLP 
Ed Depuy-Synthese , Pennsylvania , USA ) , 钢板 尺寸 
22-nim x 14 mm, 厚度 2 mm; 螺钉 直径 4 mm, J& E 
15 tm ,通过 模型 测量 并 结合 手术 实际 操作 ,本 研究 
建 这 了 椭圆 盘 型 和 箱 型 两 种 融合 器 的 初始 设计 域 结 
yz 6 所 示 。 

© 


17 mm 
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chin 


14 mm 13 mm 


15 mm 


15 mm,13 mm,6.5 mm 


(b) 箱 形 设计 域 


17 mm.14 mm.6.5 mm 


(a) 椭圆 盘 形 设计 域 


A6 两 种 融合 器 的 初始 设计 域 
Fig.6 Two kinds of primary design domains 


for the interbody fusion cage 
图 6 中 ,椭圆 盘 型 的 长 轴 为 17 mm, 短 轴 为 
14 mm, 高 为 6.5 mm; 箱 型 结构 长 宽 分 别 为 15 mm 和 
13 mm ,高 为 6.5 mm。 椭圆 盘 型 融合 器 的 生物 力学 
性 能 是 被 广泛 认可 的 ,引入 箱 型 结构 是 为 了 通过 
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比较 ,检验 本 人 研究 的 优化 结果 是 否 会 受 不 同 初始 设 
计 域 的 影响 , 即 本 拓扑 优化 方法 是 否 具 有 普 适 性 。 
将 融合 器 、 钢 板 固 定 系 统 与 上 文 验证 过 的 C3 ~ C7 
倾 椎 有 限 元 模型 进行 装配 (图 7)。 其 中 钢板 与 上 下 
椎 体 为 绑 定 连接 , 同时 融合 器 与 上 下 椎 体 的 接触 面 
也 设 为 Tie 的 连接 方式 。PEEK 材料 的 杨 氏 模 量 为 
3 660 MPa, 泊 松 比 为 0.41。 


(a) 盘 形 融合 器 植 入 模型 


(b) 箱 形 融合 器 植 入 模型 


图 7 ”两 种 构 型 融合 器 植 人 的 有 限 元 模型 


Fig.7 Finite element models for the cage insertion 


in two kinds of configurations 

1.2.2. 拓扑 优化 参数 设置 

拓扑 优化 模式 是 在 特定 的 载荷 工 况 下 , 求 体 积 
分 数 约束 下 的 最 大 刚度 设计 。 具 体操 作 在 Optistruct 
( Optistruct 12. 0, Altair) 软件 完成 。 步 又 如 下 。 

1) 创建 拓 扑 优化 设计 变量 :根据 融合 器 的 网 格 
大 小 ,实施 最 大 成 员 尺寸 约束 与 最 小 成 员 尺 寸 约 束 ， 
避免 出 现 大 面积 连续 区 域 和 小 孔洞 。 本 人 研究 划分 融 
合 器 的 网 格 大 小 为 0.5 mm, 最 小 成 员 尺 寸 设 为 
1 mm, 最 大 成 员 尺 寸 设 为 3.5 mm( 分 别 是 2 倍 和 7 
fb) 。 

2) 利 用 volumefrac 模块 ,分 别 建立 体积 分 数 和 
柔 度 响应 。 

3) 创 建 体积 分 数 约束 、 应 力 约 束 :设置 体积 分 数 
上 限 分 别 为 30% 和 50% 。 这 是 因为 文献 分 析 颈 椎 
融合 器 与 上 下 终 板 的 接触 面积 不 得 低 于 终 板 面积 的 
30% — 4096 7 ,假如 融合 器 初始 设计 域 面积 已 基本 

盖 整 个 终 板 的 支撑 区 域 , 则 优化 后 保留 30% 体积 

(接触 面 肯定 大 于 30% ) 是 可 接受 的 ,保留 50% 体积 
为 另 一 种 计算 结果 ,用 于 对 比 研究 。 

4) 定 义 目 标 函 数 :选择 最 小 化 结构 的 柔 度 (应 
变量 ) 为 目标 函数 ,使 结构 的 刚度 最 大 。 

5) 提交 Optistruct 软件 进行 优化 计算 。 
1.2.3 拓扑 优化 方法 

首先 进行 单 工 况 和 多 工 况 作用 下 的 融合 器 拓扑 
优化 ,比较 椭圆 盘 型 和 箱 型 两 种 不 同 原始 设计 域 , 选 
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择 更 优 的 结构 在 SolidWorks 软件 中 进行 二 次 优化 。 
13 拓扑 优化 融合 器 的 模拟 生物 力学 评价 


1.3.1 ACDF 手术 节 段 融合 模型 的 建立 

选取 经 SolidWorks 软件 构建 出 的 颈椎 融合 器 优 
化 模型 ,利用 布尔 运算 的 方法 建立 内 部 植 骨 块 模型 ， 
模拟 手术 中 从 骼 骨 取得 的 植 骨 材料 , 放 入 融合 器 内 
并 压 实 。 对 于 植 骨 块 的 模拟 ,本 研究 采用 三 维 实体 
单元 ,其 材料 属性 的 定义 比 松 质 骨 稍 低 , 弹 性 模 量 设 
为 100 MPa , 泊 松 比 为 0.2。 

对 融合 器 和 植 骨 模型 划分 网 格 ,网 格 单元 密度 
大 小 为 0.2, 单 个 融合 器 单元 规模 在 5 万 个 单元 左 
右 。 选 取 前 文 经 过 验证 的 颈椎 C3 ~ C7 完整 节 段 的 
有 限 元 模型 ,进行 模拟 ACDF (C5/6) 手术 的 操作 。 
其 中 ,融合 器 内 部 与 植 骨 块 进行 绑 定 约束 ,模拟 植 骨 
; 演 的 情况 。 各 融合 器 上 下 面 与 椎 骨 的 接触 面 也 为 
ARZT ,模拟 完 全 融合 的 情况 。 
1 边界 条 件 与 载荷 

" 上述 融合 模型 所 施加 的 边界 条 件 与 载荷 均 与 前 
中 完整 模型 的 验证 情况 相同 。 这 样 可 以 控制 变 
其 ( 获 完整 模型 的 计算 数据 作为 对 照 组 ,便于 与 优化 
机 曾 的 仿真 计算 数据 作 比 较 。 
e 


2>= 结 果 
>< 
LE 单 工 况 下 融合 器 优化 结果 


-S11) 前 届 工 况 下 融合 器 优化 结果 
如 图 8 所 示 。 其 中 体积 分 数 约束 为 5096 ( 简 述 

为 0.5) ;体积 分 数 约束 为 30% ( 简 述 为 0.3)。 
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图 8 两 种 构 型 融合 右 在 前 届 工 况 下 的 优化 结果 


Fig.8 The optimization results of the 


two fusion cages in flexion 
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最 终 的 优化 结果 以 单元 密度 的 形式 显示 。 红 色 
表示 单元 密度 接近 于 1 ,是 计算 后 需要 保留 的 部 分 ; 
蓝 色 表示 单元 密度 接近 于 0, 计算 后 单元 消失 ; 而 处 
于 0~1 之 间 的 中 间 密 度 单元 ,设置 单元 密度 闵 值 为 
0.5, 大 于 单元 密度 贱 值 的 部 分 保留 ,小 于 阐 值 的 部 
分 则 去 除 。 最 终 图 8 中 呈现 的 结构 实际 是 单元 密度 
BIECK T 0.5 的 部 分 。 

在 前 届 工 况 下 ,图 8(a) 表 示 当 体积 约束 分 数 为 
0.5 时 ,两 种 构 型 融合 器 的 优化 结果 ,图 8(b) 为 体积 
约束 分 数 为 0.3 时 ,两 种 构 型 融合 天 的 优化 结果 。 
总 体 优化 结果 为 舍弃 前 方 部 分 ,保留 后 方 结构 。 
为 PEEK 材料 融合 器 相对 于 颈椎 前 路 钛 板 的 弹性 模 
量 小 得 多 ,因此 前 届时 主要 由 前 路 固定 钛 板 来 承受 
载荷 , 故 融合 器 前 方 部 分 由 于 承 力 小 ,在 优化 时 单元 
自动 消失 。 

2) 后 伸 工 况 下 融合 器 优化 结果 

在 后 伸 载 荷 作 用 下 (力矩 治世 轴 负 方向 ) ,设置 
0.5 的 单元 密度 阅 值 ,得 到 图 9 的 优化 结果 。 后 伸 


时 承 力 部 位 主要 是 在 后 方 ,因此 4 个 模型 优化 结果 
都 是 保留 了 后 部 结构 ,PEEK 聚合 物 材 料 呈 层 片 状 。 


a" 
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(a) 0.5 (b) 0.3 


到 9 两 种 构 型 融合 器 在 后 伸 工 况 下 的 优化 结果 


Fig.9 The optimization results of the two 


fusion cages in extension 

3) 左 侧 及 右 侧 弯曲 工 况 下 融合 器 优化 结果 

在 颈椎 向 左 侧 和 右 侧 弯 曲 的 载荷 工 况 下 CZJAR. 
i& Z 轴 方 向 ) ,模型 优化 结果 如 图 10 ~ 图 11 所 示 。 
当 向 左 侧 弯 曲 时 ,融合 器 模型 在 0.5、0.3 体积 分 数 
控制 下 ,总 体 优化 趋势 是 在 融合 器 右 侧 形成 前 后 连 
通 的 中 空 结 构 。 当 材料 体积 分 数 增加 时 (由 0.3 增 
加 到 0. 5) ,材料 的 分 布 更 倾向 于 增加 融合 器 上 下 面 
的 表面 积 ,有 助 于 增加 与 上 下 椎 体 的 接触 面 。 


当 向 右 侧 弯 曲 时 ,其 优化 后 的 结构 与 相对 应 模 
型 在 左 侧 弯 曲 工 沉 下 的 优化 结果 基本 对 称 (图 11)。 
从 理论 上 分 析 , 人体 贷 椎 模型 沿 矢 状 面 左 右 对 称 , 如 
果 融 合 絮 的 植 入 位 置 在 正中 位 置 , 则 左 、 右 弯曲 工 况 
下 融合 器 的 受 力也 是 对 称 的 ,因此 这 两 种 工 况 下 融 
合 器 结构 优化 结果 具有 对 称 性 是 符合 逻辑 的 , 故 用 
于 模型 优化 计算 也 是 可 靠 的 。 
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Fig. 10 The optimization results of the two 
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图 11 两 种 构 型 融合 器 在 右 侧 弯曲 工 况 下 的 优化 结 


Fig.11 The optimization results of the two fusion 


cages in right bending 

4) 左 侧 及 右 侧 扭转 工 涡 下 融合 器 优化 结 

在 贷 椎 左 侧 和 右 侧 扭 转载 荷 工 况 下 (力矩 沿 Y 
轴 方 向 ) 的 优化 结果 显示 单元 分 布 集中 于 中 后 部 及 
两 侧 , 主要 是 受 前 路 钛 板 影响 。 左 右 两 侧 扭转 优化 
后 的 结构 基本 左右 对 称 ,结果 符合 逻辑 。 


2.2 多 工 况 作用 下 的 融合 器 拓扑 优化 结构 
由 于 人 体 颈 椎 的 基本 运动 分 布 在 矢 状 面 届 伸 、 
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冠状 面 侧 弯 和 水 平面 扭转 这 6 个 方向 上 ,因此 必须 
考虑 多 工 况 下 融合 融 应 满足 的 拓扑 结构 ,并 不 能 将 
每 个 模型 在 6 种 工 况 下 的 6 个 拓扑 结构 做 简单 的 县 
加 ,这 样 会 出 现 材料 元 余 , 失 去 了 优化 意义 。 所 以 本 
研究 提取 每 种 工 况 下 拓扑 结构 的 特征 进行 分 析 , 同 
时 用 加 权 和 柔 度 这 一 响应 作为 优化 目标 ,得 出 多 工 况 
下 融合 器 的 拓扑 优化 结构 。 加 权 和 柔 度 就 是 把 每 个 工 
况 下 的 柔 度 响 应 做 一 个 加 权 和 。 以 这 样 加 权 和 柔 度 为 
全 局 响应 的 拓扑 优化 结果 。 

1)6 种 工 况 下 椭圆 盘 形 融合 器 加 权 柔 度 优 化 结 
果 (PEEK 材料 ,控制 体积 分 数 不 超 过 50% 及 30% ) 
如 图 12 所 示 。 

图 12 是 同一 融合 屁 在 不 同体 积分 数 约束 下 的 
两 个 优化 结果 。 第 一 行 是 优化 后 保留 的 结构 ,其 单 
元 密度 在 0.5 及 以 上 ,下 一 行 图 中 颜色 接近 蓝 绿色 
的 部 分 是 单元 密度 小 于 0.5 的 部 分 ,在 优化 后 单元 
消失 ,将 其 提取 出 来 放 在 图 中 便于 观察 。 图 12 显 
示 : 当 体积 分 数 为 0.3 时 ,会 使 前 方 结构 缺失 ,材料 
集中 分 布 于 后 部 及 两 侧 ; 在 体积 分 数 分 别 为 0.5 和 
0.3 时 ,融合 咒 优 化 为 多 孔径 和 单一 大 孔径 的 形状 ; 

前 者 比 后 者 增加 了 表面 积 。 
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Fig.12 Topology structure of elliptic cage 
2)6 种 工 况 下 PEEK 材料 箱 形 融合 器 加 权 柔 度 
优化 结果 ( 控制 体积 分 数 不 超 过 0.5 及 0.3) :与 前 
文 椭圆 盘 型 融合 器 的 优化 结果 一 致 (图 13)。 
根据 图 12 和 图 13 中 的 优化 结果 ,用 OSSmooth 
模块 根据 单元 密度 结果 生成 封闭 的 等 值 曲 面 ,这 样 
得 到 拓扑 优化 构 型 的 几何 表面 ,可 用 于 进一步 的 
CAD 建 模 。 图 14 中 模型 类 别 采 用 “体积 分 数 约束 ” 
的 方式 来 表示 。 
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(c) 箱 形 0.5 
图 14 融合 器 的 几何 表面 


Fig. 14 Geometric surface of the fusion cages 


(d) 箱 形 0.3 
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2.32 融合 器 二 次 优化 与 优化 结果 提取 


“三 前 文 已 提 及 ,引入 椭圆 盘 型 和 箱 型 两 种 初始 设 


na 
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融合 器 最 终 的 几何 模型 ,如 图 15 所 示 。 


(a) 0.5 (b) 0.3 
图 15 两 种 盘 形 融 合 器 设计 构 型 


Fig.15 Two configuration designs of the elliptic cage 


2.4 拓扑 优化 融合 器 的 生物 力学 评价 结果 


2.4.1 颈椎 各 节 段 活动 度 
在 预 载荷 与 1.5N * m 力矩 作用 下 ,有 贷 椎 在 3 个 

运动 面 上 产生 运动 。 各 个 节 段 在 力矩 作用 下 运动 到 
一 个 平衡 位 置 ,这 个 最 终 的 平衡 位 置 距离 初始 位 置 
转动 了 一 定 角度 ,把 相 邻 两 个 节 段 间 的 旋转 角度 之 
差 称 为 关节 活动 度 ,也 即 相 对 旋转 角度 。 经 计算 ,得 
到 各 个 节 段 在 矢 状 面 的 前 届 / 后 伸 的 关 广 活动 度 之 
和 ,冠状 面 左右 侧 弯 的 活动 度 之 和 ,以 及 水 平面 左右 
扭转 的 活动 度 之 和 ,数据 如 表 2 所 示 。 

表 2 优化 模型 在 屈伸 、 侧 弯 、 扭 转载 荷 下 的 

各 节 上 段 活 动 度 ( 单位 :*) 
Tab.2 The range of motion of optimized model 


in flexion-extension , bending and torsion (unit:°) 


Ht, 379. REIS FR BALL 8 Fd d UC MR 
否 会 有 不 同 的 变化 趋势 。 而 上 述 结果 已 表明 ,本 研 
究 采 用 的 拓扑 优化 技术 并 不 受 构 型 变化 的 影响 ,两 
种 初始 设计 域 表现 出 了 相同 的 优化 结果 。 尽 管 早期 
的 生物 力学 研究 发 现 箱 型 融合 器 相 比 其 他 形态 的 融 
合 器 ,能 够 在 贷 椎 各 个 方向 的 运动 中 , 均 提供 足够 的 
稳定 性 “” ,但 是 由 于 箱 型 融合 器 尖锐 的 8 个 棱角 有 
损伤 相 邻 椎 体 终 板 的 风险 , 故 理想 的 融合 器 四 周 应 为 
圆 弧 形 ,这 与 目前 市 场 上 的 融合 器 构 型 基本 相符 。 

由 于 前 文 已 经 判定 融合 器 初次 优化 后 一 定 会 出 
现 大 的 中 空 结构 , 故 在 二 次 优化 中 ,首先 在 融合 器 内 
挖 除 一 个 直径 为 7 mm 的 中 央 孔 ,将 这 一 融合 器 模型 
FLAME C5 ~ C6 节 段 ,用 与 上 文 相同 的 拓扑 优化 
方法 对 其 进行 结构 拓扑 优化 。 直 接 采 用 加 权 和 柔 度 优 
化 方法 ,最 终 在 SolidWorks 软件 中 建 出 了 椭圆 盘 型 


项 目 C3 ~ C4 C4 ~ C5 C5 ~ C6 C6 ~ C7 
前 届 + 后 人 
正常 模型 7.30 8.63 7.9 7.61 
Peek-0. 5 6.73 7.94 0.26 7.55 
Peek-0. 3 6.73 7.94 0.26 7.55 
双 侧 弯曲 
正常 模型 6.47 6.17 4. 64 4.93 
Peek-0. 5 6.26 6.00 0.18 4.89 
Peek-0. 3 6.26 6. 00 0.18 4.89 
双 侧 旋转 
正常 模型 7.82 6.86 5.10 4.61 
Peek-0. 5 7.70 6.65 0.19 4.60 
Peek-0. 3 7.70 6.65 0. 19 4. 60 


由 表 2 可 知 , 手 术 后 融合 节 段 C5 ~ C6 的 关节 活 
动 度 显著 降低 ,不管 是 届 伸 、 侧 弯 还 是 扭转 ,都 降低 到 
0.3" 以 下 ,相对 于 正常 模型 4" ~ 8° 的 活动 度 ,数据 降 
低 了 90% 以 上 。 在 不 同 构 型 的 融合 器 模型 间 , 非 手术 
节 段 活动 度 差异 很 小 ,与 正常 模型 相 比 略 有 降低 。 


2.4.2 融合 器 与 植 骨 块 的 应 力 分 布 

Chiang 等 报道 当 颈 椎 届 曲 时 , 由 于 前 路 钛 板 
会 分 担 更 多 的 应 力 ,对 融合 需 形 成 应 力 遮 挡 , 会 降低 
其 内 部 植 骨 块 所 受 的 应 力 。 而 当 颈 椎 后 伸 时 ,前 路 
铁 板 的 固定 作用 会 削弱 ,使 融合 需 面 临 更 大 的 应 力 ， 
其 至 有 发 生 移 位 和 下 沉 的 风险 。 故 本 研究 集中 对 不 
同 构 型 融合 器 在 后 伸 运 动 中 的 应 力 分 布 结果 进行 分 
析 ,如 图 16 所 示 。 


(b) 9E Zr 450.3 


(c) 植 骨 块 0.5 


(d) 植 骨 块 0.3 
图 16 融合 器 与 植 骨 块 的 应 力 分 布 图 
Fig. 16 The stress on the cage and bone graft 
图 16 表明 :优化 后 的 融合 器 整体 应 力 值 较 小 ， 
其 按 缘 的 应 力 高 于 内 部 ;而 内 部 植 骨 材 料 的 承受 应 
加 起 对 较 大 ,后 部 的 应 力 高 于 前 部 。 
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日 本 研究 对 不 同 初始 构 型 不同 体积 分数 控制 的 
融合 絮 模 型 进行 拓扑 优化 分 析 , 目标 是 在 特定 的 生 
理 载 荷 下 ,寻求 有 限 材 料 的 较 优 分 布 , 即 在 确保 融合 
器 自身 稳定 性 (结构 刚度 ) 的 前 提 下 ,增加 其 柔 度 ， 
以 便 有 效 减少 应 力 遮 挡 , 使 融合 器 内 部 的 植 骨 块 能 
够 分 担 更 多 的 应 力 ,促进 融合 。 在 这 样 的 拓扑 方案 
设计 下 ,本 研究 根据 拓扑 优化 后 得 到 的 结果 ,提出 了 
2 种 融合 器 材料 构 型 。 上 述 2 种 构 型 的 融合 器 结构 
有 一 个 共同 特点 , 那 就 是 内 部 中 空 ,材料 主要 集中 在 
外 围 区 域 。 这 样 的 构 型 不 仅 适 合 植 骨 ,也 符合 贷 椎 
节 段 间 的 力学 环境 。Grant 等 ” 的 实验 表明 , 终 板 
后 部 与 外 部 所 承担 的 载荷 明显 强 于 前 部 与 中 部 ， 
此 合理 的 融合 器 支撑 区 域 并 不 在 中 央 , 而 是 在 外 于 
和 后 部 , 故 本 研究 中 的 融合 器 设计 是 合理 的 。 

椭圆 盘 型 和 箱 型 2 种 构 型 的 融合 器 ,表现 出 了 
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相同 的 优化 结果 。 按 形状 类 别 基 本 可 分 为 两 类 : 单 
一 大 孔径 和 多 孔径 融合 器 。 当 材料 体积 分 数 为 
3096 I] , 模型 为 基本 的 环形 ; 当 材 料 体积 分 数 为 
5096 时 ,放弃 单一 大 孔径 的 设计 ,在 内 部 中 空 的 基础 
上 增加 融合 器 上 下 面 的 表面 积 ,并 通过 小 和 孔 与 中 空 
部 分 连通 ,便于 植 骨 。 这 种 筛 状 的 构 型 是 本 研究 优 
化 的 特征 结果 ,能 有 效 增加 融合 器 与 植 骨 床 的 接触 
面积 ,同时 又 不 影响 移植 骨 与 植 骨 床 的 接触 生长 。 
本 研究 中 拓扑 优化 的 融合 器 会 显著 增加 与 相 邻 终 板 
的 接触 面积 ,有 助 于 防止 出 现 融 合 器 下 沉 。 

计算 机 模拟 的 生物 力学 结果 显示 C576 节 段 的 
活动 基本 被 限制 ,是 因为 前 路 铁 板 的 固定 作用 以 及 
优化 后 的 融合 器 仍 可 提供 足够 的 支撑 强度 ,使 得 相 
应 节 段 难以 产生 相对 运动 。 因 此 ,不 同 构 型 融合 器 
在 活动 度 方面 并 没有 显著 的 差异 ,配合 前 路 猴 板 国 
定 的 效果 相当 。 融 合 器 所 受 应 力 的 下 降 是 去 除 材 料 
宛 余 和 加 权 柔 度 优化 的 结果 。 此 消 披 长 ,这 也 造就 
了 植 骨 块 内 部 应 力 的 增加 ,有 助 于 促进 融合 器 与 相 
邻 终 板 的 融合 。 


4 结 论 


对 于 适用 于 ACDF 手术 的 PEEK 材料 颈椎 融合 
带 , 经 过 有 限 元 辅助 下 的 三 维 拓扑 优化 后 ,结构 理 
想 ,生物 力学 性 能 满意 。 结 论 如 下 。 

1) 本 研究 的 优化 方法 具有 普 适 性 ,即使 对 于 不 
同 的 初始 构 型 ,也 能 达到 相同 的 优化 结果 。 

2) 在 一 定 的 体积 分 数 约束 下 ,融合 器 整体 结构 为 
环形 ,通过 进一步 加 权 和 柔 度 优化 使 其 前 部 结构 减少 ， 
甚至 缺失 ,材料 集中 分 布 于 融合 器 的 后 部 及 两 侧 。 

3) 得 状 构 型 的 融合 器 能 有 效 增 大 融合 器 与 相 邻 
终 板 的 接触 面积 ,理论 上 可 以 降低 终 板 所 受 应 力 ,能 
否 有 助 于 防止 融合 器 下 沉 , 还 有 待 于 更 多 的 生物 力 
学 试验 证 实 。 

4) 优化 后 的 颈椎 融合 器 通过 合理 去 除 宛 余 材料 
和 增加 有 柔 度 , 既 能 扩大 植 骨 容量 ,又 能 有 效 减少 应 力 
遮挡 ,利于 植 骨 块 的 应 力 传导 ,将 有 助 于 促进 融合 。 
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